A multicolor near- infra red study of the dwarf nova IP Peg 
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^ , ABSTRACT 

We report the analysis of JHK S light curves of the eclipsing dwarf nova IP Peg in quiescence. The light 
curves are dominated by the ellipsoidal variation of the mass-donor star, with additional contributions from 
the accretion disc and anisotropic emission from the bright spot. A secondary eclipse is visible in the J and 
H light curves, with 2% and 3% of the flux disappearing at minimum light, respectively. We modeled the 
observed ellipsoidal variation of the secondary star (including possible illumination effects on its inner face) 
to find a mass ratio of q = 0.42 and an inclination of i = 84.5°, consistent in the three bands within the 
uncertainties. Illumination effects are negligible. The secondary is responsible for 83%, 84% and 88% of the 
flux in J, H and K s , respectively. We fitted a black body spectrum to the JHK S fluxes of the secondary star 
to find a temperature of T b b = 3100 ± 500 K and a distance of d — 115 ± 30 pc to the system. We subtracted 
the contribution of the secondary star and applied 3-D eclipse mapping techniques to the resulting light curves 
to map the surface brightness of a disc with half-opening angle a and a circular rim at the radius of the bright 
spot. The eclipse maps show enhanced emission along the stream trajectory ahead of the bright spot position, 
providing evidence of gas stream overflow. The inferred radial brightness-temperature distribution in the disc 
is flat for R < 0.3Rli with temperatures ~ 3500-ftT and colors consistent with those of cool opaque radiators. 
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1. Introduction IP Peg is a relatively bright (V ~ 14.5 mag), long period 
(f) ■ (Porb = 3.8 hr) dwarf nova showing recurrent outbursts ev- 
. ■ The shape of the secondary star in cataclysmic variables ery few months in which the system increases its brightness 
(CVs) and low-mass x-ray binaries (LMXBs) is defined by by „ 2 mag m the visible Tne binary is seen almost edge on 
^ . its Roche equipotential surface. Modeling the ellipsoidal ( inc i matj ion i > 80°) allowing for eclipses of both the white 
variation produced in the light curves of these binaries by dwar f/ acC retion disc and the secondary star, as well as of 
the changing aspect of the distorted secondary star with bi- the bright spot (hereafter BS) . Optical and ultraviolet (UV) 
nary phase provides constraints on the orbital parameters light curves of IP Peg are dominated by anisotropic emis- 
(e.g., Warner 1995) and yields an estimate of the contribu- sion from the prorninent B S, making it hard to constrain 
tion of the secondary star to the total light. the system parame ters (e.g., Wood & Crawford 1986) and, 
Unlike visible and ultraviolet light curves, which are therefore, difficult to apply eclipse mapping techniques to 
dominated by emission from the disc and bright spot, the investigate its accretion disc structure. Szkody & Mateo 
near-infrared (NIR) light curves of CVs and LMXBs have (1986) found evidence of ellipsoidal variations in JHK S 
significant contributions from the cool secondary star. The light curves of IP Peg. Froning et al. (1999) [hereafter F99] 
determination of the secondary star flux in a set of wave- f oun d a pronounced double-hump modulation after sub- 
lengths allows a direct estimate of its spectral type and trading the ellipsoidal variation of the secondary star from 
also opens up the possibility of inferring the distance to their H band light curve. Their eclipse mapping analysis 
the binary by means of photometric parallax. Furthermore, indicates that the accretion disc is cool with a flat radial 
knowing the contribution of the secondary star to the light temperature distribution (T bb ~ 3000 K) in the H band, 
curve allows one to isolate the light from the accretion disc 

and to apply eclipse mapping techniques (e.g., Baptista & This paper reports a multicolor study of IP Peg in the 
Steiner 1993) to investigate the disc structure in the in- JHK S bands. Section [2] presents the observations and data 
frared. reduction. Section [3] describes the procedure used to fit the 
ellipsoidal variation and the application of 3-D eclipse map- 
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in memoriam ping techniques to the light curves after we subtract the 



contribution of the secondary star. The results are discussed 
in Section |4] and summarized in Sectional 



2. Observations and data reduction 

IP Peg was observed with WHIRCAM (Hughes et al. 1996) 
at the 4.2m William Herschel Telescope in La Palma on 
1996 October 26-29, while the system was in quiescence. 
The observations comprise 3 orbital cycles in the H band, 
and a bit less than 1 orbital cycle each in the J and K s 
bands. All runs were performed in good weather with bright 
moon, but no clouds. The seeing ranged from 1.0" to 1.6". 
The exposure times were of l-3s (H) and of 3s (J and K s ). 
The observations are summarized in Table [TJ 

We used a dithering procedure to estimate the contri- 
bution of the sky background, nodding the telescope in a 
five-position square pattern (center plus four corners) . The 
sky level was obtained from the median of the five images 
of each set and was then subtracted from each image of the 
set. A series of dark images, with the same exposure time 
as the science frames, were obtained during the run at in- 
tervals of ~ 1 — 2 hrs. These images were used to remove 
the dark current of the chip. Correction of flat-field effects 
was also performed. 

Data reduction was performed using APPHOT/IRAF0 
aperture photometry routines. The frames of each dithering 
sequence were aligned with an interactive procedure based 
on the IRAF 'register' task and combined to increase the 
signal-to-noise ratio (S/N) of the measured stars. The non- 
linearity of the detector was corrected using the 'irlincor' 
routine in the CTIO package. Fluxes were then extracted 
for the variable and for a comparison star 7.4" North and 
7.5" East of the variable. 

Differential light curves (target star flux divided by com- 
parison star flux) were computed. The light curves were flux 
calibrated using the 2MASS JHK S zero point constants 
(Skrutskie et al. 2006) and the absolute magnitudes for the 
comparison star (J = 9.82, H — 9.56 and K s = 9.48). 
No color term corrections were applied. The data were 
phase-folded according to a modified version of the linear 
ephemeris of Wolf et al. (1993), 

T mid {HJD) = 2445615.4156(4) + 0.15820616(4) • E, (1) 

were T m id is the time of inferior conjunction of the sec- 
ondary star. The resulting light curves can be seen in Fig. 
[TJ The H band data were combined to improve the S/N of 
the light curve. The J and K s band data do not cover the 
full orbit of the binary. 



3. Data analysis 

3.1. The white dwarf eclipse width A</>. 

The optical study of Wood & Crawford (1986) leads to a 
white dwarf eclipse width in the range A</> = 0.086 — 0.092. 
The uncertainty in estimating A(j) arises from the fact that 
the white dwarf eclipse ingress in IP Peg is veiled by the 
much more pronounced ingress of the BS eclipse and by the 



Table 1. Journal of the observations. 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 
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Fig. 1. Light curves of IP Peg with modeled ellipsoidal vari- 
ation added to the inferred (constant) disc contribution to 
fit the data. Vertical dashed lines mark the ingress/egress 
phases of the white dwarf for an eclipse width of A(j> = 
0.0918 cycles. 



large amplitude flickering seen in optical light curves prior 
to BS ingress. 

In order to refine the value of Ac/), we employed an iter- 
ative procedure, assuming binary parameters derived from 
the fit of the ellipsoidal variation (Sect. I3-2|) and applying 
the corresponding phase offset <po needed to make the ob- 
served white dwarf mid-egress feature coincident with phase 
Acj)/2. The procedure searches for the pair of (A</>, </) ) val- 
ues that yields the best % 2 for the model light curve. This 
iterative process converges to a solution with A(j> = 0.0918 
and 4>o = +0.01, consistent with the upper limit of Wood 
& Crawford (1986). We use the J band light curve in this 
process because it shows clearly the white dwarf egress, un- 
like the H band light curve. The K s band light curve was 
not used for this procedure because of its lower S/N and 
reduced phase coverage. Because the white dwarf egress 
feature is present in both J band eclipses, we may safely 
exclude the possibility that the observed feature is due to a 
flicker/flare from the disc. The light curves shown in Fig.Q] 
were corrected by the derived value of <po in order to make 
mid-eclipse coincident with phase zero. 



3.2. Fitting the ellipsoidal variation. 



Table 2. Modeled IP Peg parameters 



The ellipsoidal variation is modeled with the aid of a light- 
curve synthesis program, which computes the flux emitted 
by a Roche-lobe filling star as its aspect changes with binary 
phase. The surface of the Roche lobe filling star is divided 
into a large number of tiles. The flux emitted by each tile is 
modified to account for gravity- and limb-darkening effects. 
We used the gravity-darkening coefficient of Sarna (1989), 
(3 = 0.05, and adopted the non-linear square-root limb- 
darkening law of Diaz-Cordoves & Gimenez (1992), 

I = Io(l-a(l-coay)-b(l-y/mFr)), (2) 

where 7 is the angle between the line of sight and the normal 
to the surface of the tile, and Iq is the intensity emitted 
normal to the tile surface. We assumed aj = —0.465, an = 
-0.454, a K = -0.448 and bj = 1.199, b H = 1.173, b K = 
1.066 (Claret 1998). 

In order to account for both limb and gravity darkening 
we need to adopt values for the temperature T po i e and grav- 
ity g po ie at the pole of the star. Using a procedure similar 
to that described by F99, we tested a range of physically 
plausible values varying T po i e from 2800 to 3200^ and logg 
from 3.5 to 5.0. The results were rather insensitive to the 
choice of parameters, except for the K s band in which the 
program only converges for values of T po i e = 3000 — 3200X 
and \ogg — 4.5. We then decided to adopt T po i e — 3000-ftT, 
log g — 4.5 and solar metallicity, since these values are in 
better agreement with those in the literature (Szkody & 
Mateo,1986; Leggett, 1992; F99). 

The parameters fitted by the ellipsoidal variation code 
are: mass ratio (q), binary inclination (i), the peak intensity 
of a spot at the inner hemisphere of the secondary star 
relative to the mean intensity at the star's surface (I s /Iq), 
the orientation angle of the spot with respect to the line 
joining both stars (p s ), a phase-independent additive flux 
level , and the flux of the secondary star at phase zero (F s ). 
The program minimizes the \ 2 between the observed light 
curve and the ellipsoidal variation model for a given set of 
parameters with an amoeba minimization scheme (Press et 
al. 1986). Alternatively, it is possible to fit a phase offset 
that minimizes the x 2 of each set of parameters (see Sec. 

ECU- 

The NIR light curves of IP Peg show contributions from 
other sources apart from the ellipsoidal variation of the sec- 
ondary star. As a first step, we removed only the phases 
covering the primary and secondary eclipses from the light 
curve before attempting to fit the ellipsoidal variation. This 
leads to hard-to-converge, unrealistic solutions which un- 
derestimate the hump at phase +0.75 and overestimate the 
hump at phase +0.25 in the J and H bands. The difference 
in brightness is larger in the J band - where the asymmetry 
in the eclipse shape caused by the BS is larger - and disap- 
pears in the K s band, which shows almost no evidence of 
the BS in the eclipse shape. This indicates that anisotropic 
emission from the BS contributes significantly to the hump 
centered at phase +0.75. We therefore decided to omit the 
phase range [+0.6; +0.9], as well as the primary and sec- 
ondary eclipses, from the fitting procedure. With the new 
restriction, it was possible to fit a model light curve to the 
data of all bands. 

Table [2] lists the results of the fitting procedure. The 
model light curves are shown as solid lines in Fig.[TJ Our re- 
sults indicate that the secondary star is responsible for 84% 
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F s is the flux of the secondary and Ft is the total flux 



and 88% of the total brightness in the H and K s bands, re- 
spectively, in agreement with the results of Szkody & Mateo 
(1986) and F99. Littlefair et al. (2001) fitted an M4V type 
star to the K s band spectrum of IP Peg to find that the sec- 
ondary star contributes 62% of the total light at that wave- 
length. However, as pointed out by Harrison et al. (2005a, 
b), the depth of the absorption lines from the secondary star 
in CVs are reduced with respect to isolated stars of similar 
spectral type, perhaps by some atmospheric effect. In this 
case, attempts to match the depth of the absorption lines 
will systematically underestimate the contribution of the 
secondary star to the total light. This effect may account 
for the lower contribution inferred by Littlefair et al. (2001). 
Alternatively, the different inferred contributions could be 
a consequence of changes is disc brightness with time, with 
a slightly brighter disc leading to a lower secondary star 
relative contribution at the epoch of the observations of 
Littlefair et al. (2001). It is hard to test this possibility 
because IP Peg is at the limit of detection for amateur 
astronomer while in quiescence, showing a typical scatter 
of 0.2 mags in its historical light curve. We analyzed the 
AAVSO historical light curve of IP Peg and found no dis- 
cernible (larger than 0.2 mag) difference in brightness state 
between the epochs of the observations of F99, Littlefair et 
al (2001) and this work. 

The fitted Gaussian spot on the inner face of the sec- 
ondary star is roughly centered at the LI point and gives a 
negligible contribution to the total flux in all bands within 
the uncertainties, indicating that irradiation effects are not 
significant in the IR continuum. 

An inclination versus mass ratio diagram for IP Peg 
is shown in Fig. [21 where the derived range of values of 
q and i for the three bands are depicted. There is good 
agreement between the results from the three bands, with 
the dispersion of the mass ratio values being much smaller 
than the formal errors of the determinations in the J and 
H bands. The best fit (lower x 2 value) is obtained for the H 
band data. Because this light curve is the average of data 
from three orbits, it has the most complete phase cover- 
age, higher S/N, and lower influence from flickering. The 
J and K s light curves include data from only one orbit 
with incomplete phase coverage, although the lack of phase 
coverage of the J band data does not affect the fit of its 
ellipsoidal variation. In spite of their smaller formal error, 
the results for the K s band are less reliable since only a 
small part of the light curve was used for the fit. Taking 
into account the constraint derived from the inferred width 
of the white dwarf eclipse (Sect. 13. ip . our best solution is 
i = 84.5° and q = 0.42. This set of parameters is indicated 
by a filled circle in Fig. fJl 
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Fig. 2. Inclination - mass ratio diagram. Dashed lines in- 
dicate the relations for A(f> = 0.0863 (upper curve) and 
A</> = 0.0918 (lower curve) and bracket the range of (i,q) 
values of Wood & Crawford (1986). Vertical dotted lines 
mark the range of values determined by Watson et al. 
(2003). Three boxes indicate the 1 — a limits on the de- 
rived parameters from the JHK S band data. A cross, an 
open circle, an open diamond and a filled circle depict the 
formal solution of Wood & Crawford (1986), Beekman et 
al (2000), Marsh (1988) and this work, respectively. 



Figure [2] also compares our results with those in the 
literature. The binary parameters derived from the ellip- 
soidal variations in the J and H band data are in rea- 
sonably good agreement with those of Wood & Crawford 
(1986), Beekman et al. (2000) and Watson et al. (2003). 
The solution of Marsh (1988) relies on inferring the white 
dwarf radial velocity from the wings of the emission lines 
from the disc. This technique is prone to large uncertainties, 
particularly when there is significant phase offset between 
the spectroscopic inferior conjunction of the secondary star 
and the observed mid-eclipse time - as it is in the case in 
IP Peg. 

Stellar atmosphere models (e.g. Kurucz 1979) and avail- 
able stellar atlases (e.g. Strecker et al. 1979) are not suitable 
for fitting IR fluxes of very cool stars, either because of a 
poor spectral resolution and sampling in the IR or because 
the range of temperatures of the grid/atlas does not include 
stars cooler than ~ 3300X. Therefore, we fitted black bod- 
ies to the extracted JHK S fluxes of the secondary star of IP 
Peg in order to estimate its temperature and distance. The 
best fit solution yields Tbb — 3100 ± 500-FT and a distance 
of dbb = (115 ± 30)pc, assuming R2 = 0.4i? Q (Beekman et 
al. 2000). The results are consistent with those of Szkody 
& Mateo (1986). 

3.3. 3-D eclipse mapping 

After subtracting the ellipsoidal variation corresponding 
to the adopted binary parameters from each light curve, 
eclipse mapping techniques (Baptista & Steiner 1993) were 
applied to the residual curves to derive maps of the sur- 



face brightness distribution of the IP Peg accretion disc m 
JHK S . 

Our eclipse map is a 3D surface consisting of a 51 x 
51 pixels disc grid with a half-opening angle a (the angle 
between the midplane and the disc surface) covering the 
primary Roche lobe (up to R — Rli) plus a circular rim 
of 101 pixels at R = Rb s , the inferred radius of the BS. 
With this geometry there is no need to extract the orbital 
hump from the light curve (F99), since it can be accounted 
for by emission from the circular rim. We also note that 
this geometry allows the disc brightness distribution to ex- 
tend beyond the radial position of the BS. The radius of 
the BS was derived as follows. We measured the BS mid- 
ingress/egress phases (4>bi = —0.027, (t>b e = +0.983) in the 
J band light curve after correcting for 4>o (Sect. 13. lj) . For 
the adopted binary geometry (i,q), the pair of (4>u,4>be) 
values maps into an x-y position in the orbital plane which 
consistently falls along the ballistic stream trajectory. The 
BS radius is taken as the radius of the circle that passes 
through this position, Rb s = 0.58Rli- The adopted binary 
geometry and the BS radius are depicted in Fig. 01 For the 
eclipse mapping modeling, only the data in the phase range 
[-0.15;0.15] was analyzed. 

We performed Monte Carlo simulations to estimate the 
uncertainties of the reconstructions (Rutten et al. 1992). 
For a given light curve a set of 100 artificial light curves 
is generated, in which the data points are independently 
and randomly varied according to a Gaussian distribution 
with standard deviation equal to the uncertainty at that 
point. The light curves are fitted with the eclipse mapping 
code to produce a set of randomized eclipse maps. These 
are combined to produce an average map and a map of 
the residuals with respect to the average, which yields the 
statistical uncertainty at each pixel. The uncertainties ob- 
tained with this procedure are used to draw the contour 
maps of Fig. and to estimate the uncertainties in the de- 
rived radial brightness-temperature distributions (Fig. [5]) 
and the flux-ratio diagram (Fig. [6]) . 

The disc half-opening angle a is a free parameter in the 
problem. The entropy of the eclipse map is a useful tool 
in gauging the correct value of a. Overestimating (under- 
estimating) the disc half-opening angle introduces spurious 
structures in the disc side closest to (farther from) the LI 
point. Because the entropy is a measurement of the smooth- 
ness of the eclipse map, these structures are flagged with 
higher entropy. Therefore, one may estimate the value of a 
by performing a set of reconstructions for a plausible range 
of a and selecting the one with lowest entropy. Simulations 
(Borges, Baptista & Catalan 2007) show that it is indeed 
possible to use the entropy as a criterion to infer the value 
of a, but that the map of lowest entropy underestimates 
the correct a by 1.5 — 2.0°, depending on the binary geom- 
etry {i,q). We confirmed this finding with careful simula- 
tions done with the specific geometry of IP Peg, for which 
we find the offset to be 1.5°. In order to infer the value 
of a we used the K s band light curve, which has the most 
symmetric eclipse shape and the smallest contribution from 
the orbital hump. We tested a range of a values between 
0° — 8° with increasing steps of Aa = 0.5°. The map of low- 
est entropy and highest degree of symmetry is obtained for 
a = 0.5°. According to the above simulations, we adopted 
a = 2° for the eclipse mapping reconstructions. This is 
consistent with computations of vertical disc structure by 
Meyer & Meyer-Hofmeister (1982), Smak (1992) and Hure 




Fig. 3. Left: The light curves after correction of ellipsoidal variation (dots) and the best fit 3-D eclipse mapping model 
(solid curve). Vertical dotted lines mark the ingress/egress phases of the white dwarf. Middle-left: Corresponding disc 
brightness distributions on a logarithmic gray-scale. Dotted lines show the Roche lobe. A solid contour line is overploted 
on each eclipse map to indicate the 3-cr confidence level on the intensities. The logarithmic intensity scale is indicated by 
the horizontal gray-scale bar (in erg/cm 2 /s/Hz/str). Middle-right: Asymmetric disc brightness distributions. Dotted 
lines show the Roche lobe, the ballistic stream trajectory and a disc of radius 0.58i?Li. Right: Disc rim brightness 
distributions. The x-axis indicates the azimuth with respect to the line joining both stars. A vertical dotted line marks 
the azimuthal position of the bright spot. The dotted lines show the X-a limits on the intensities. 



& Galliano (2000), which predict a ~ 1.5° - 2.0° for mass 
accretion rates M = 10~ n - 10~ 10 M Q yr _1 . 

Light curves, corresponding eclipse maps and disc rim 
brightness distributions are shown in Fig. [31 For a better 
visualization of structures in the disc brightness distribu- 
tions the asymmetric disc components are also shown. A 
symmetric component is obtained by slicing the disc into a 
set of radial bins and fitting a smooth spline function to the 
mean of the lower half of the intensities in each bin. The 
spline-fitted intensity in each annular section is taken as the 
symmetric disc-emission component. This procedure pre- 
serves the baseline of the radial distribution while removing 
all azimuthal structure. The asymmetric disc component is 
then obtained by subtracting the symmetric disc from the 
original eclipse map (Saito & Baptista 2006). The asym- 
metric disc component accounts for 30%, 20% and 16% of 
the total flux in J, H and K s , respectively. 

The maps show enhanced emission along the ballistic 
stream trajectory close to the white dwarf, providing ev- 



idence for gas stream overflow in IP Peg. The maximum 
of the enhanced emission is wavelength-dependent, occur- 
ring farther downstream in J than in K s , possibly because 
of the progressively larger amount of gravitational energy 
released in collisions as the stream approaches the white 
dwarf. This tail of enhanced emission along the ballistic 
stream ahead of the BS position is reminiscent of that seen 
in the dwarf nova WZ Sge (Spruit & Rutten 1998; Skidmore 
et al. 2000). Spruit & Rutten (1998) point out that such a 
tail is to be expected as a consequence of the post-impact 
hydrodynamics of the stream. 

A BS is expected to form at the intersection of the bal- 
listic stream with the disc outer edge. The ballistic stream 
hits the 0.58i?Li disc rim at an azimuth 8 = 15.6° with 
respect to the line joining both stars 0. The JHK S disc rim 
distributions (Fig. [3j show a bright spot at an azimuth of 

2 The difference between the value of 8 for the adopted q=0.42 
geometry and for the alternative q=0.5 case (Wood & Crawford 
1986) is negligible. 
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Fig. 4. The geometry of IP Peg. A filled circle marks the 
expected position of the BS; an open circle depicts the ob- 
served azimuth of maximum emission along the disc rim. 
An arrow indicates the direction of maximum emission of 
the orbital hump, tilted forward 54° from the azimuthal 
position of the bright spot. The large circle corresponds to 
a disc of radius 0.58i?Li- 



9 ~ 20° with azimuthal extent (full width half maximum) 
A8 ~ 30°, in reasonable agreement with the predicted po- 
sition of the BS. The intensity of the spot decreases with 
increasing wavelength, in accordance with the correspond- 
ing reduction in the height of the orbital hump. There is 
also a tendency for the centroid of the spot to move towards 
shorter azimuths with increasing wavelengths. 

It is interesting to compare the azimuthal position of 
the BS in the disc rim map with that of the orbital hump 
maximum. The azimuth of hump maximum is 6 ~ 74° 
(<j) m ax — —0.2, Fig. [l|), or ~ 54° forward from the azimuth 
of the BS. Figure |4] shows a schematic diagram of IP Peg, 
where the offset between the azimuth of the BS and that of 
the maximum hump emission is clear. A similar effect has 
been seen in other dwarf novae. If BS emission is produced 
in a shock between disc and stream gas, the maximum emis- 
sion will be normal to a direction between the disc and 
stream flows (Warner 1995, p. 81). For example, in OY Car 
the hump maximum is displaced ~ 19° forward with respect 
to the azimuth of the BS (Wood et al 1989). The offset in 
IP Peg is larger. Here, the azimuth of orbital hump maxi- 
mum coincides with the normal to the gas stream at the BS 
position, suggesting that the impact shock lays along the 
stream trajectory with negligible influence from disc ma- 
terial. We note that the phase of maximum of the orbital 
hump in the NIR-light curves does not coincide with the 
one in the optical (cj> — 0.9, see Fig. 1 of Wood & Crawford 
1986), as already noticed by F99. 

Figure \E\ shows the JHK S disc radial brightness tem- 
perature distributions for an assumed distance of d = 115 
pc. Dotted lines show the T cx i? -3 / 4 law of opaque steady- 



state discs for mass accretion rates of M= 10 , 10 and 
10 -12 MQj/r -1 , assuming Mi — 1.05M Q (Beekman et al. 
2000) and Ri = 0.0124i?© (derived using the white dwarf 
mass-radius relation of Nauenberg 1972). 

The inner disc (R < O.SRli) shows a flat brightness 
temperature distribution reminiscent of those seen in qui- 
escent dwarf novae (e.g., Wood et al. 1986, 1989), with in- 
ferred temperatures of ~ 3500 K at R = O.IRli for the 
JHK S bands within the uncertaities. The consistency of 
the inferred brightness temperatures in the three bands 
suggests that the inner disc is optically thick with emis- 
sion close to blackbody. The temperatures decrease in the 
outer disc regions (R > 0.3Rli) in reasonable agreement 
with the T cx i?~ 3 / 4 law in the J and K s bands and at a 
lower shallow gradient in the H band. There is marginal ev- 
idence of higher brightness temperatures in the outer disc 
regions in the H band with respect to the J and K s bands. 
If real, such H band excess would indicate that the gas in 
the outer disc regions is opaque with a vertical temperature 
gradient (such that the H~ free-free and bound-free opac- 
ity minimum at 1.6 microns leads to a relative increase in 
outcoming flux in the H band with respect to the J and 
K s bands). This is hard to reconcile with the results of 
Littlefair et al (2001), which suggest that the outer regions 
of the accretion disc in IP Peg are optically thin, unless the 
source of the observed mirror eclipses is an optically thin 
chromosphere above the (optically thick) accretion disc. 

A flux-ratio diagram of the accretion disc of IP Peg is 
shown in Fig. [6j The colors for the symmetric component of 
the eclipse maps are plotted together with relationships for 
blackbody (BB), main-sequence stars (MS) and optically 
thin HI emission (HI). The BB and HI spectra were com- 
puted with the synphot/IRAF package. The BB and HI 
fluxes are extracted by convolving the respective spectrum 
with the response function of each infrared passband, and 
the flux ratios are computed. The MS colors were extracted 
from Bessell & Brett (1988) and transformed to the 2MASS 
photometric system using the relations of Carpenter (2001). 
The uncertainties in the disc colors are quite large and in- 
crease towards the outer (and fainter) disc regions. The un- 
certainties in the colors are dominated by the lower S /N of 
the K s band light curve and eclipse map. The colors of the 
symmetric inner disc (R<0.3 Rli) are consistent with those 
of cool opaque radiators (with T co i or ~ 3000 — 5000 K) at 
the 1-cr confidence level. The disc becomes redder (cooler) 
with increasing radius and starts to deviate from the BB re- 
lationship for R > 0.3i?Li. The observed trend in outer disc 
colors is consistent with that inferred from the brightness 
temperature distributions and again suggests an iJ-band 
excess flux for the outer disc regions. The disc colors move 
toward the upper left of the diagram, in the direction op- 
posite to that expected for optically thin gas. While the 
error bars are large and the diagram should be viewed with 
caution, the results suggest that the outer disc of IP Peg is 
opaque with a vertical temperature gradient. 

3.4. The secondary eclipse 

Figs. Q] and [7| shows the model curve from the ellipsoidal 
variation together with the data. It is clear from the J and 
H-band data that some occulting medium is covering part 
of the light from the secondary star between phases 0.45 

3 The observed difference in temperature is at the lcr limit. 
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Fig. 5. Left: Disc temperature brightness-distribution on 
a log scale for the symmetric disc component. The solid 
curve connecting the filled squares is the radial profile 
for the disc and the dashed curves indicate the 1-cr lim- 
its obtained with Monte Carlo simulations. Dotted lines 
show opaque steady-state disc models for mass accretion 
rates of M= 1CT 10 , lCT 11 and 10~ 12 M e yr~ l , assuming 
Mi = 1.05 M and Ri = O.O124R . Right: Disc rim 
brightness temperature distribution on a log scale. A dotted 
line indicates the azimuthal position of the bright spot. 



and 0.55. Eclipses of the secondary star by the disc were 
seen in the if -band data of F99 - who pointed out that 
the shallow depth of the eclipse requires that part of the 
intervening accretion disc must be optically thin - and also 
in the spectroscopic data of Beekman et al. (2000) - who 
deduced that the accretion disc of IP Peg must be opaque to 
account for the depth of the secondary eclipse. Littlefair et 
al. (2001) found evidence of an eclipse of the secondary star 
by optically thin parts of the accretion disc from their K- 
band time-resolved spectroscopy of IP Peg in quiescence. 
Assuming local thermal equilibrium, they found that the 
occulting optically thin gas is quite hot (10000 K < T < 
20000 K) and accounts for, at least, the outermost 20 per 
cent of the disc. Here we model the secondary eclipse to 
infer the colors of the occulted secondary star face and to 
estimate the radial extent of the opaque occulting disc. 

Assuming the binary geometry derived in section 13.21 
we simulated the eclipse of the secondary star by an opaque 
cylinder of radius rd and half-opening angle a and searched 
for the values of r d and a which minimize the % 2 of the 
fit to the data. In doing this we are making the simplistic 
assumption that the disc is fully opaque up to a certain 
radius and becomes fully transparent thereafter. If opacity 
changes with radius and height above the disc mid-plane 
one expects that a sizeable fraction of the disc will partially 
transmit the light from the secondary star and a larger 



region would be needed to occult the same amount of light. 
Therefore, our exercise provides a rough lower limit to the 
size of the opaque occulting disc. Having this in mind, the 
best-fit solutions are r d (J) = 0.19±0.06 R L1 , a(J) = (0.7 ± 
0.4)° (corresponding to an occultation of 2.3 ± 0.4 per cent 
of the secondary star light) and r d (H) — 0.22 ± 0.03 Rli, 
a(H) = (1.5 ±0.7)° (occultation of 3.0 ±0.6 per cent of the 
secondary star light), respectively, for the J and H bands. 
The secondary eclipse and the best-fit model curves are 
shown in Fig. [71 The model provides a good description of 
the data for phases <f> < 0.55 but deviates after the eclipse 
because of the (unaccounted) contribution of the orbital 
hump to the light curve. Assuming a disc radius of 0.58 
Rli, this suggests that at least the inner ~ 1/3 of the 
accretion disc is opaque in IP Peg in quiescence. 

We fitted a blackbody model to the extracted J and H 
band deficit fluxes at mid-secondary eclipse to find a tem- 
perature of the occulted, inner hemisphere of the secondary 
star of T oc = 2300 ±600 K. After correcting for the gravity 
darkening effect, this becomes T' oc = 2500 ± 700 K, consis- 
tent at the la limit with the blackbody temperature of the 
(outer hemisphere of the) secondary star derived in section 
13.21 The agreement between the inferred temperatures of 
the inner and outer hemispheres of the secondary star is in 
line with the conclusion derived from the ellipsoidal varia- 
tion fit, namely, that illumination effects on the secondary 
star of IP Peg in quiescence are negligible in the NIR con- 
tinuum. The occulting area corresponds to > 13 per cent 
of the projected surface of the secondary star at phase 0.5. 

4. Discussion 

F99 found a double-humped modulation in their i/-band 
light curve after subtraction of the ellipsoidal variation from 
the secondary star (see their Fig. 5). There is no evidence 
of double-humped modulation in our data. Since the phas- 
ing of their data was not secured by a clear identification 
of the white dwarf egress feature, their analysis may be af- 
fected by errors in the assumed binary phases. F99 added a 
phase offset of = +0.027 to their light curves to correct 
for deviations of the ephemeris. However, precise measure- 
ment of white dwarf eclipse egress times from contemporary 
observations of IP Peg in quiescence indicate a phase off- 
set of A(j) = +0.008 for the epoch of the F99 observations 
(Baptista et al. 2005). Thus, the light curves of F99 are sys- 
tematically shifted by 8<j> = +0.019 towards positive phases. 
Indeed, there is an upward kink in the flux at phase ~ +0.06 
in their data (best seen in their Fig. 11) which seemingly cor- 
responds to the white dwarf egress, suggesting a phase offset 
~ +0.015 with respect to the expected white dwarf egress 
phase. Simulations show that such a small phase offset be- 
tween the data and the ellipsoidal variation model would 
be enough to introduce a spurious double-humped modu- 
lation in the light curve after removal of the contribution 
from the secondary star. The phase offset also displaces the 
whole disc brightness distribution towards the trailing side 
of the disc (the one containing the gas stream) and would 
account for much of the offset between the expected posi- 
tion of the BS and the position derived from their eclipse 
map. 

The modeled ellipsoidal variation also indicates that the 
irradiation of the secondary star surface is negligible in the 
infrared. This result is in apparent contradiction with those 
of Davey & Smith (1992) and Watson et al (2003). Despite 
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Fig. 6. Flux-ratio diagram. A solid line indicate the colors 
of black body emitters. Illustrative temperatures along the 
sequence are indicated by labels. A dashed line shows the 
colors of optically thin emitting gas with three tempera- 
tures labeled. A dotted line indicates stellar main-sequence 
colors, with some spectral types labeled. Filled squares con- 
nected by solid lines indicate the color of the symmetric disc 
component. The error bars represent the standard deviation 
with respect to the mean value in each radial bin. Labels 
mark the radius in units of Rli- 



the differences in the results, both works find a significant 
decrease in Nal A8190A line strength over the inner face of 
the secondary star in IP Peg as a consequence of irradiation 
effects. We remark that, if the irradiated energy does not 
penetrate deeply in the atmosphere but mainly heats the 
upper atmospheric layers, it may reduce the vertical tem- 
perature gradient (and, therefore, lead to a decrease in Nal 
absorption line strength) without affecting the continuum 
radiation arising from deeper atmospheric layers. This is 
in line with the investigation of irradiation effects on CV 
secondaries by Barman, Hauschildt & Allard (2002). Their 
preliminary study shows that, for a typical CV, irradiation 
leads to a significant change in the temperature structure 
of the uppermost atmospheric layers of the secondary star 
(where the Nal line comes from) leaving the deeper layers 
(the r = 1 region where the IR continuum is produced) 
mostly unaffected. 

Our disc temperatures are higher than those inferred 
by F99. We find temperatures of ~ 3500if in all bands 
while they find temperatures of ~ 3000iC in H . This is 
in line with the fact that our 77-band light curve has 2.1 
mJy more flux than the 1994 Sep light curve of F99 (at 
the reference phase 4> — 0.25, where the contribution of 
BS anisotropic emission should be minimal). Our data and 
those of F99 were collected, respectively, 5 and 6 weeks after 
an outburst. The observed difference in disc temperatures 
could be explained if the accretion disc of IP Peg slowly 
cools down during quiescence. Alternatively, the observed 
difference in light curve flux and disc temperatures could 
be accounted for by the 10% uncertainty in flux calibration 
of both works. 
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Fig. 7. J and -ff-band light curves of IP Peg (open sym- 
bols) with modeled ellipsoidal variation without (dashed 
line) and with (solid lines) secondary eclipse included. 



5. Summary 

We modeled the ellipsoidal variations in JHK S light curves 
of IP Peg to estimate the system parameters and the contri- 
bution of the secondary to the total light (Table [2|) . Using 
the calibrated fluxes we inferred an equivalent blackbody 
temperature of 3100 ± 500 K, and we derived a distance to 
the system of 115 ± 30 pc. 

Using 3D eclipse mapping techniques we derived the 
surface brightness distribution of the IP Peg accretion disc. 
With the aid of the entropy of the eclipse map, we find 
a disc half-opening angle of a = 2°. The disc brightness 
distribution shows an asymmetric structure along the gas 
stream trajectory in the inner disc regions close to the white 
dwarf in all bands (R ~ 0.1i?£i), indicating the existence 
of gas stream overflow in IP Peg. 

The amplitude of the orbital hump decreases with in- 
creasing wavelength. The hump is modeled as an extended 
(A(9 ~ 30°) bright spot located at the edge of the accretion 
disc with brightness temperatures of ~ 10000X in J and 
H and ~ 6000-K" in the K 8 band. The phase of maximum 
of the orbital hump in NIR-light curves, <j> = 0.74, does not 
coincide with the one in the optical (</> = 0.9, see Fig. 1 
of Wood & Crawford 1986) nor does it coincide with the 
phase of maximum emission of a radially emitting bright 
spot. The eclipse position of the bright spot differs from its 
theoretical azimuthal position by ~ 5° and its direction of 



maximum emission (inferred from the light curve) differs 
by ~ 54° from the direction of radial emission. 

In order to match the depth and width of the secondary 
eclipse, a sizeable fraction of the accretion disc must be 
opaque. The inferred temperatures of the inner and outer 
hemispheres of the secondary star are the same to within 
the uncertainties, indicating that illumination effects are 
negligible for the IR continuum. 
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